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Considering a variety of dynamical scenarios within the lin- 
ear a model, we examine the conditions for pionic modes to 
be amplified as the chiral field relaxes towards the normal vac- 
uum following the transient restoration of approximate chiral 
symmetry in a high-energy collision. While Bjorken-type lon- 
gitudinal expansions appear to be insufficient for amplifica- 
tion for occur, analoguos expansions in two or three dimen- 
sions may enhance the low-energy power spectrum by up to 
an order of magnitude. 

PACS numbers: 25.75.+r,11.30.Rd,11.30.Qc,12.38.Mh 



The possibility of forming disoriented chiral conden- 
sates (DCC) in high-energy collisions has generated con- 
siderable research activity in recent years (for a recent 
review, see rcf. 0). The basic premise is that such colli- 
sions produce extended hot regions within which approx- 
imate chiral symmetry is temporarily restored. The sub- 
sequent non-equilibrium relaxation may then lead to the 
formation of coherent sources of low-energy pions and as- 
sociated anomalous multiplicity distributions ^ ||. Im- 
portant insights into the dynamics of DCC formation 
have been gained with analytical approaches |5]-|J and, 
in the most elaborate studies, the evolution of the chi- 
ral degrees of freedom has been simulated numerically 
|5] Q. We present here a simple and instructive frame- 
work within which we analyze the conditions for the oc- 
currence of the DCC phenomenon. 

Most dynamical studies of disoriented chiral conden- 
sates have been based on the linear a model in which the 
chiral degrees of freedom are described by the real 0(4) 
field (p(r, t) = (a, 7r) having the equation of motion 



[□ + \(4> 2 - v 2 )] <t> = Hcr . 



(1) 



The three parameters in the model can be fixed by spec- 
ifying the pion decay constant, /„■ = 92 MeV, and the 
meson masses, — 138 MeV and m a — 600 MeV, 
leading to the values A = (ml - m%)/2f* = 20.14, 
v = [(ml - 3ml)/(ml - m^)] 1 / 2 ^ = 86.71 MeV, and 
H = mlU = (120.55 MeV) 3 , with K,c=l @. 

As is apparent from eq. (|l|) , the vacuum configuration 
is aligned with the a direction, </> vac = (/„., 0), and at low 
temperature the fluctuations represent nearly free a and 
7T mesons. In the other extreme, at temperatures well 
above v, the field fluctuations are centered near zero and 
approximate 0(4) symmetry prevails. 

In order to analyze the situation, it is instructive to 
decompose the chiral field, (f>(r, t) — 4>(t) + S(p(r, t). The 
first term, <fi, is the average over a suitable region in space 



(of dimensions larger than the correlation length) and 
can be identified with the (local) order parameter, while 
the fluctuations, 5<j>(r), may be considered as elementary 
quasi-particle excitations. We may simplify the discus- 
sion without affecting the conclusions by assuming that 
the order parameter is always fully aligned, <f> = (<tq, 0). 
The fluctuations are then 6</>(r,t) = (8a, Sir). 

Applying a Hartree-type factorization [jl5 17 1, we may 
replace the full equations of motion (jl]) by a set of self- 
consistent approximate equations, 



[□ + Mo] ct = H , 
[□ + /4] 5a = , 
[□ + /4] Sir = , 

where the associated effective masses are given by 
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2 < 5a 2 > 
2 < Sa 2 > 
2 < Sn 2 > 



(2) 
(3) 
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with Sni being the particular cartesian component of Sir. 
Thus the field fluctuations provide an additional stiffncs 
resisting the growth of the order parameter ao. The effec- 
tive masses are degenerate for <7o=0 and they vanish at 
the temperature T c = y/2v, in the mean-field treatment. 
Moreover, we always have jUq < /j, 2 < fi 2 . 

Eqs. @ and (||)were first derived in ref. JTsj . The term 
< 5cj> 2 > in eqs. (Jqfl) 1S ^ ne sum °f ^ ne ne ld fluctuations 
in each of the N=A chiral directions and is of leading or- 
der in 1/N. These are the 'direct' terms that have been 
included in previous DCC treatments or discussions in 
terms of effective masses [E|B Jl(]i[RO | . The next term 
in eqs. (^0) is twice the fluctuation along the particu- 
lar direction considered (either parallel or perpendicular 
to the order parameter) and arises from the 'exchange' 
terms. A recent analysis of the statistical properties of 
the linear a model ]lj suggests that their effect is sig- 
nificant (~ 2/V=50%) and that their inclusion leads to 
a quite good approximation to the full eq. ([j]). 

The approximate equations (||-Q) provide a convenient 
framework for developing a qualitative understanding of 
the dynamics generated by the full equation ([!]). Imag- 
ine that the system is initially created in thermal equi- 
librium at a temperature To well above T c . The field 
fluctuations are then sufficiently large to ensure fi 2 >0 in 
all three eqs. (|^-0). The system is expected to experience 
a cooling resulting from expansion and radiation, so the 
fluctuations decrease in the course of time. This reduces 
f/r which allows the order parameter to grow larger, thus 
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counteracting the decrease of the effective masses. The 
resulting behavior of /i 2 is then a delicate balance: for 
slow cooling the induced growth of Co is approximately 
adiabatic and the system relaxes towards the vacuum 
through metastable configurations; however, if the fluc- 
tuations diminish rapidly a compensating growth of the 
order parameter can no longer occur quickly enough and 
one or more of the effective masses may turn imaginary, 
/x 2 <0, indicating that the system has entered a regime 
exhibiting exponential growth of some modes. 

Fig. |l| shows the equilibria and the unstable region. 
As T is increased, the fluctuations grow steadily and the 
equilibrium value of <7o decreases from / w . The most 
rapid change occurs at T «220 MeV, above which <7rj 
tends to zero. The effective pion mass \i v increases mono- 
tonically with T from its free value towards « 1.6T 
for r>T Cl while \i a first decreases, then displays a min- 
imum at T w 240 MeV, and finally becomes degenerate 
with fi n . The border of the unstable region intersects the 
a axis at /„■ and extends up to T c /v3 at <to=0. 
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T = 400 MeV 



D=0: sphere 
D=1 : matter 
D=2: matter 
D=3: matter 
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FIG. 1. Dynamical trajectories. 
The combined dynamical evolution of the order param- 
eter er =< & > and the field fluctuations < S<j} 2 > 1 / 2 
(sampled for r <2.5 fm in the case of D=0). The dashed 
curve connects the equilibria from T—0 to above 500 MeV 
and the unstable region within which /i 2 < is shown 
by the shaded region. Each system has been prepared in 
thermal equilibrium at To = 400 MeV, using a periodic 
box (20 fm side length). The irregular solid trajectory 
(-D=0) was obtained by solving the standard eq. (|l|) af- 
ter applying a spherical Saxon- Woods modulation factor 
(5 fm radius and 0.5 fm width) to the hot matter, thereby 
producing a hot sphere embedded in vacuum jl{J. The 
other three trajectories have been obtained without a 
spatial modulation but with the term —(D/t)dt4> added 
in the equation of motion to emulate uniform expansions 
in D=l,2,3 dimensions. The marks along the trajectories 
are positioned at time intervals of At = 0.2 fm/c. 



Also shown in fig. [t] is the dynamical trajectory of the 
central part of a Ni-sized spherical source prepared at 
T =400 MeV. The system keeps away from the unstable 
regime, exhibiting an approximately adiabatic evolution. 
(The windings reflect the damped oscillations of <7q in 
the relaxing effective potential, in the familiar manner 
pyr^f .) This behavior is rather robust, as it occurs for 
a wide range of initial temperatures and for rod or slab 
geometries as well. It thus appears that initially static 
field configurations in local equilibrium do not develop 
any instabilities during their subsequent expansion. 

However, it is expected that the early parton dynamics 
causes the chiral field to be in a state of rapid expansion. 
The subsequent evolution may then lead to a supercooled 
configuration situated inside the unstable region, thus ef- 
fectively producing a "quench" . A number of quenched 
scenarios have been considered but they were 

imposed by fiat, thereby reducing the predictive power 
of the dynamcial calculations (essentially any degree of 
magnification can be achieved by suitable adjustment of 
the initial conditions). It is our aim here to reduce the 
degree of arbitrariness by elucidating under which condi- 
tions a quench-like early scenario may develop dynami- 
cally from plausible initial configurations. 

Simple Bjorken-like pictures can be invoked to emu- 
late expanding scenarios in D dimensions, either longi- 
tudinal (D=l) |n]|n^,[l(|, transverse (D=2), or isotropic 
(D=3) jl2|,^0|]. We have considered such scenarios in the 
following approximate manner. The essential effect of 
the expansion is the appearance of an additional term 
on the right-hand side of eq. (|l|), —(D/t)d t (f>, where the 
time variable should now be reinterpreted as the elapsed 
proper time in a comoving frame (starting at to — 1 fm/c, 
usually). The corresponding Lorentz transformation of 
the (scaled) spatial coordinates is less essential for our 
present discussion and has therefore been ignored. 
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FIG. 2. Idealized isotropic expansions. 
Similar to fig. |l], but for idealized isotropic expansions 
using D='i in the modified equation of motion. 
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The additional term in the equation of motion, 
— (D/t)dt<p, acts as a time-dependent damping, its form 
being akin to the Rayleigh dissipation function in clas- 
sical mechanics, and it is significantly more effective in 
reducing the fluctuations than self-generated expansions. 
In order to examine its effect, we ignore the spatial geom- 
etry and consider a macroscopically uniform configura- 
tion within a large box. The discussion is then simplified 
and the resulting scenarios can be regarded as idealized 
representations of chiral matter subjected to an exter- 
nally prescribed cooling rate, and so the results will have 
a corresponding general applicability. 

Figure |l] includes dynamical trajectories obtained in 
this manner for 13=1,2,3. The effect increases with D, 
since the dimensionality of the expansion effectively acts 
as the strength of the damping term. The isotropic ex- 
pansion leads to a significant incursion into the unstable 
region, while the longitudinal one is too slow for that. 
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Elapsed time t (fm/c) 
FIG. 3. Time evolution of the pion mass. 
The time evolution of the square of the effective 
pion mass, eq. (Q), for idealized D-dimensional expan- 
sions starting from thermal equilibrium at To=400 MeV. 
The thin curves show the corresponding evolution of /!q 
(recall < /zj < fi^r). All curves approach the free pion 
mass TOjr (dotted line) in the course of time. 



TABLE I. Amplification coefficients/correlation lengths. 



To (MeV) 


D = 1 


D — 2 


D = 3 


200 


0.00 /1.4 


0.02 /I.8 


0.11 /2.0 


220 


0.00 /1.3 


0.50 /1.9 


0.55 /2.5 


240 


0.01 /1.3 


1.20 /2.0 


1.19 /2.7 


300 


0.00 /0.9 


1.84 /1 .7 


2.06 /2.7 


400 


0.00 /0.6 


1.67 /1.3 


2.49 /2.1 


500 


0.00 /0.5 


I.3I /I.I 


2.61 /1.6 



The enhancement factor G£ =0 in macroscopically uni- 
form chiral matter and the resulting correlation length, 
expressed as the full width at half maximum (in fm) of 
the spatial two-point function -< 5tv(x) ■ 5ir(y) >~, for 
pions emerging after idealized expansions with £>=1,2,3 
starting from thermal equilibrium at Tq. 



Figure H shows trajectories for Z?=3 starting from vari- 
ous temperatures. If the initial temperature is lower than 
200 MeV or so, the initial value of do is already fairly 
large (over 60 MeV) and the dynamical trajectories will 
miss the unstable region. A wide range of higher temper- 
atures lead into the unstable region, provided the super- 
cooling occurs sufficiently fast. Ultimately, at very high 
temperatures (above those shown) the system will again 
stay stable throughout, because it takes so long to re- 
duce the fluctuations down to critical size that the order 
parameter has meanwhile had time to start its growth. 

In order to quantify the analysis, it is useful to consider 
the time evolution of the effective masses. Since \j? a > /ij 
we concentrate on fi* which is illustrated in fig. |[ It 
is noteworthy that /ip w /ij throughout the evolution. 
(Early on <ro is small so 0(4) symmetry holds approxi- 
mately, and later on the fluctuations are less important 
than Co-) The amplification of the lowest pionic modes is 
then practically identical to that of the order parameter 
itself. This simple feature makes it an easier task to ana- 
lyze more complicated scenarios. As was already evident 
from figs. 1-2, multiple incursions into the unstable re- 
gion are possible, especially for large initial cooling rates 
(since a rapid quench brings <7o into larger oscillations 
around f^), but the first one is always dominant. 

It is convenient to express the resulting enhancement 
of a mode in terms of its amplification coefficient plj , 



G7T 
Ic 



dt y/-Uk(t) 3 



(8) 



where the frequency is given by the dispersion relation, 



= k 2 



t4, 



using the time-dependent effective pion 



mass (Q). [The dispersion relations for the quasiparticles 
are in principle affected by the presence of the expan- 
sion term, as is ordinarily the case for a damped oscil- 
lator. This effect may be significant for damping rates 
large enough to bring the system into the unstable re- 
gion, but has been ignored in the present study] The 
quantity exp(GJ) expresses approximately the factor by 
which the amplitude of pions having the wave number k 
has been magnified due to the incursion(s) into the un- 
stable regime. Since largest magnification occurs for k=0 
(a finite k adds a positive amount to w|), we concentrate 
on this quantity as an upper bound. Moreover, we note 
that the minimum in fJ-^(t) (see fig. |3|) determines the 
maximum wave number for which magnification occurs. 
Amplification coefficients obtained for various idealized 
expansion scenarios are shown in table ||, together with 
the resulting width of the pion correlation function. 

The purely longitudinal expansions (D=l) miss the 
unstable region, except for insignificant incursions near 
T=240 MeV, while a significant degree of magnifica- 
tion occurs for the transverse and isotropic expansions, 
amounting to over a factor of ten in the most favorable 
cases. It should be recalled that Gq applies to k—0 only, 
so it provides an upper bound on the enhancement. 
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This qualitative sensitivity to the collision dynamics 

1 ' ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ^ highlights the importance of employing realistic initial 

r =400 MeV - conditions for the dynamical simulations of DCC forma- 

Z tion. Ultimately the appropriate initial field configura- 

D ° o o d=1 I tions must be calculated on the basis of the early partonic 

° ++ + + ° + D=2 z evolution, a task which is thus crucial for our ability to 

j + , * '- )= ^ I assess the prospects of forming disoriented chiral conden- 

+ + 6» 8 ta - sates in high-energy collisions. 

OS g 
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FIG. 4. Pion power spectrum. 
The relative power spectrum of the pions, ~ ^k^k' 
where 7Tfc is the Fourier amplitude of the pion field, plot- 
ted as a function of the pion kinetic energy tOk-m^. The 
extraction is made at large times when the asymptotic 
scenario of free evolution has been reached. The plots 
are based on samples of 20 field configurations prepared 
at To =400 MeV and subjected to idealized expansions 
with either D—l (open dots), D=2 (crosses), or D=3 
(solid dots). The irregularities are primarily due to the 
shell structure in the level density of the cube. 

Figure ^ gives an impression of the net effect on the 
power spectrum of the emerging free pions. As expected, 
the transient instabilities present for D=2,3 lead to sig- 
nificant enhancements of the power carried off by soft 
pions. The effect amounts to about an order of magni- 
tude for D=3 (relative to the smooth spectrum obtained 
for D=l), in accordance with the amplification coefficient 
given in Table 0. Although these results were calculated 
for idealized expansion scenarios, they do support the 
suggestion that such enhancements may provide an ob- 
servable signal of DCC formation 22 1. 

It should be emphasized that the dynamical evolution 
of the chiral field was obtained by solving the full equa- 
tion of motion (jl|). The simple Hartree approximation 
(HI?]) has been invoked only as a convenient framework 
for the discussion since it brings out so clearly the con- 
ditions for amplification. 

Our analysis shows that the occurrence of instabili- 
ties, and the associated amplification of pionic modes, 
depends sensitively on the cooling rate, which in turn 
is intimately related to the character of the expansion. 
Our idealized scenario for D=3 corresponds closely to 
the isotropic expansion considered in r cfs. and our 

results corroborate the conclusion in |12| that such a sce- 
nario leads to amplification. Furthermore, our analysis 
suggest that a longitudinal expansion alone is insufficient 
to cause a quench, if the initial fluctuations are of ther- 
mal magnitude. This is consistent with what was found 
in refs. |11 1M for effectively one-dimensional expansions. 
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